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Abstract: Stability of linear systems with uncertain bounded time- varying delays 
is studied under assumption that the nominal delay values are not equal to zero. 
An input-output approach to stability of such systems is known to be based on the 
bound of the L2-norm of a certain integral operator. There exists a bound on this 
operator in two cases: in the case where the delay derivative is not greater than 1 
and in the case without any constraints on the delay derivative. In the present note 
we fill the gap between the two cases by deriving a tight operator bound which is an 
increasing and continuous function of the delay derivative upper bound d > 1. For 
d oo the new bound corresponds to the second case and improves the existing 
bound. As a result, delay-derivative-dependent frequency-domain and time-domain 
stability criteria are derived for systems with the delay derivative greater than 1. 
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1. INTRODUCTION 

Two main approaches have been applied to sta- 
bility analysis of linear systems with uncertain 
time- varying delay: a direct Lyapunov approach 
and an input-output approach (see e.g. Gu et al., 
2003) , which reduces the stability analysis of the 
uncertain system to the analysis of the class of 
systems with the same nominal part but with 
additional inputs and outputs. In the existing 
literature the uncertain time-varying delay has 
been divided into two types: the slowly-varying 
delay (with delay derivative less than d < I) and 
the fast-varying delay (without any constraints 
on the delay derivative) (see e.g. Kolmanovskii 
& Myshkis , 1999; Niculescu , 2001). Recently a 
third type of moderately varying delay has been 
revealed in Fridman & Shaked (2005), where the 



delay derivative is not greater than 1 (almost for 
all t). This has been obtained by applying the 
input-output approach to stability. It is known Gu 
et al. (2003), Kao & Lincoln (2004) that the latter 
approach to systems with time-varying bounded 
delays is based on the bound of the L2-norm of a 
certain integral operator. 

In the present paper we fill the gap between the 
case of the delay derivative not greater than 1 and 
the fast-varying delay by deriving a new integral 
operator bound. This bound is an increasing and 
continuous function of the delay derivative bound 
d > 1. In the limit case (where d — > c») which 
corresponds to the fast-varying delay, the new 
bound improves the existing one. As a result, 
improved frequency-domain and time-domain sta- 



bility criteria are derived for systems with the 
delay derivative greater than 1. 

Notation: Throughout the paper the superscript 
'T' stands for matrix transposition, 7?." denotes 
the n dimensional Euclidean space with vector 
norm || • ||, 7^"^™ is the set of all n x m real 
matrices, and the notation P > 0, for P e 7^"^" 
means that P is symmetric and positive definite. 
The symmetric elements of the symmetric matrix 
will be denoted by *. L2 is the space of square 
integrable functions v : [0, 00) — > C" with the 
norm \\v\\l^ = [J^ \\v{t)\\^dt]^^'^ , \\A\\ denotes the 
Euclidean norm of a n x n (real or complex) matrix 
A, which is equal to the maximum singular value 
of A. For a transfer function matrix of a stable 
system G{s), s € C 

||G||oo= sup \\G{iw)\\, t = V^. 

— OO<10<OO 



2. PROBLEM FORMULATION 

We consider the following linear system with un- 
certain time- varying delay T{t) : 

xit)=Aax{t)+Aixit~T{t)), (1) 

where x{t) G 7?." is the system state, Ai, i = 0,1 
are constant matrices. 

The uncertain delay T{t) has a form 

T{t) = h + rj{t), \vit)\<f,<h, (2) 

where /i is a known nominal delay value and /i is a 
known upper bound on the delay uncertainty. In 
the existing literature Kolmanovskii & Myshkis 
(1999), Niculescu (2001), Gu et al. (2003) the 
following types of uncertain time-varying delays 
are usually considered: 

Case A (slowly- varying delay): T{t) is a dif- 
ferentiable almost everywhere function, satisfying 



f{t)=r]{t) <d=l+p, 
where — 1 < p < 0; 



(3) 



Case B (fast- varying delay): T{t) is a measur- 
able (e.g. piecewise-continuous) function. 

Recently a moderately-varying delay with f{t) < 
d = 1 was introduced in Fridman & Shaked 
(2005). In the present note we enlarge the latter 
class of delays as follows: 

Case C (moderately- varying delay): r(i) is a 
differentiable almost everywhere function, satisfy- 
ing (3) with p > 0. 

In the present note we will improve the stability 
results in cases B and C by applying input-output 
approach and by deriving new inequalities. The 
results are easily generalized to the case of any 
finite number of the delays. 



We represent (1) in the form: 
x{t) = AQx{t) -h Aix{t -h)- 



J-h-r] 



x{t + s)ds. 

(4) 

Following Fridman & Shaked (2005) we introduce 
the following auxiliary system: 



x{t) = Aox{t) + Aix{t -h)+ iJ,Aiu{t), 
y{t) = ^/W^x{t), 



(5) 



with the feedback 



u{t) 



1 



y{t + s)ds, (6) 

where : [— l,oo] i?"*" is a scalar function 
which will be shortly defined and p is given by (3). 
The results for the delay of case B correspond to 
p = 00, i.e. to J-"(oo) in the input-output model 
(5), (6). Substitution of (6) in (5) readily leads to 
(4). 

We are looking for J^{p) which satisfies the follow- 
ing inequality 



HI 



< 



\y 



\l2' 



V 2/ e L2[0,oo), y\ 



-00,0] 



0, 
(7) 

where u is given by (6). This is equivalent to the 
fact that \x\/ Tip) is an upper bound on the Li- 
norm of the integral operator A : L2[0,oo) 
i2[0,oo) 



z{t) 



Ay{t) = |_ 



h—r} 



y{t + s)ds, y\ 



-00,0] 



i.e. that 



\l2 



\\z\\h<f^'Hp)\\y\\] 

VyeL2[0,oo), 



0. 



0, 
(8) 

(9) 



Our objective is to find T{p) (as small as possible) 
such that (7) (or equivalently (9)) holds. 

For — 1 < p < (case A) it was established in 
Gu et al. (2003) that T{p) can be chosen to be 

1. For p > the following was found in Fridman 
& Shaked (2005): T{0) = 1 and T{p) = 2 for 

(0,00]. 

We note that the value 1 of T{p) for — 1 < p < 

can not be improved (i.e. chosen to be less than 
1). Indeed, taking constant delay r]{t) = jj, which 
satisfies the condition of case A for any — 1 < p < 
0, we consider the functions ye{t) = I a.s < t < 6 , 
and yeit) = a.s t > 0. Using formula (6) with 
T{p) = 1 we immediately obtain 



\\ye\\U = o'^ Mi 

and hence ||u||L2/||j/e||L2 



1 as 9 ^ 00. 



In the present paper we will improve the values 
of T{p) for p > by showing that !F{p) can 

be chosen as a continuous increasing function of 
p > satisfying T{0) = 1 (as in Fridman & 



Shaked (2005)), but T{p) < T{oo) = 1.75 for 
p > 0. The improved values of J^{p) will readily 
lead to improved stability criteria. 



3. MAIN RESULTS 
3.1 New Bounds 

Proofs of the Lemmas of this section are given in 
the Appendix. 

Lemma 1. Consider case C. For all y{t) G 
i2[0, oo) and such that y{t) = Vt < and for 
u{t) given by (6) inequality (7) holds with given 

by 



2p+l 

p+1 ' 
fp- 1 

4p ' 



if < J9 < 1, 
if P > 1 • 



(10) 



As it was mentioned above, J-" is increasing contin- 
uous function satisfying for p > the following in- 
equality: 1 = < J^{p) < limp^oo J'ip) = 7/4. 

Lemma 2. Consider case B. For all y{t) € 
L2[0, oo) and such that y{t) = < and 
for u{t) given by (6) inequality (7) holds with 
jr(oo) := 7/4. 

Remark 1. The value 7/4 = 1.75 for JF(oo) in 
Lemma 2 is not far from an optimal one. The 
following example shows that it cannot be less 
than 1.5. Namely, define scalar functions y{t) and 
7?(i) by 



t, if < t < /i, 

— t, t£ 1^ <t < 2/i, 
0, if y{2ii -y)<0, 



vit) = 



—fj,, if t < /u, 
fj,, if t > ^. 



Setting in (6) J^(oo) = 3/2 we have u{t) = 



z{t + h)= [ y{s)ds 

Jt-ri{t) 



'-(i + M)V2, 
-(/z^ -t- 2iJLt-2t^)/2, 
= < {6fj,t-3fj,^ -2t^)/2, 
{t - 3/x) V2, 

We achieve equality in (7) since 



if - fi<t + h<0, 
ifO <t + h<iJ,, 
if /i < t + /i < 2/i, 
if 2/i < i + /i < 3/i, 
otherwise. 



2 
3/i' 



2 M 



3.2 A Tight 
Criterion 

We assume 



Frequency-Domain Stability 



Al Given the nominal value of the delay h > 0, 
the nominal system 



x{t) = Aox{t) + Aix{t - h), 
is asymptotically stable. 



(11) 



The auxiliary system (5) can be written as y = Gu 
with the transfer matrix 



G{s) = ^/W<^sl{sl - Ao - Aie-''«)-V^i- 

(12) 

By the small gain theorem (see e.g. Gu et al. 
(2003) the system (1) is input-output stable (and 
thus asymptotically stable, since the nominal sys- 
tem is time-invariant) if ||G||oo < 1. A stronger 
result may be obtained by scaling G: 

Theorem 1. Consider (1) with delay given by (2). 

Under Al the system is asymptotically stable if 
there exists non-singular matrix X such that 

||Gx||oo<l, Gx{s)=XG{s)X-\ (13) 

where G is given by (12) with T{p) of (10) and 
where p G [0,00) corresponds to case C, while 
.F(oo) = 7/4 corresponds to case B. 

Remark 2. From Theorem 1 it follows that under 
Al (1) is asymptotically stable if 

1 , . 1 



\sI{sI-Ao -Aie-'*^)-Ui| 



By Fridman & Shaked (2005) T{p) = 2,p > 
and thus (1) (with T{t) < 1 > or with T{t) 

of case B) is asymptotically stable for r(t) £ [h — 
/i, h+fj] , where /i < 0.7071fc. By the new boimds of 
Lemma 2 and Lemma 1 we obtain a wider stability 
intervals: 

p = 0.1, f{t) < 1.1, J^{p) = 1.0909, n < 0.9574fc, 

p = 0.5, f {t) < 1.5, T{p) = 1.3333, /i < 0.8660fc, 
P=l, f {t) < 2, J^Ip) = 1.5, II < 0.8165fc, 
p = 00, case B, J^{p) = 1.75, /x < 0.7559fc. 

(14) 



3.3 On Improved Stability Criteria in the 
Time-Domain 

By applying the time-domain results of Fridman 
& Shaked (2005) via descriptor model transfor- 
mation with the corresponding simple Lyapunov- 
Krasovskii functional we obtain: 

Theorem 2. System (1) is asymptotically stable 
for all delays of (2), if there exist nxn matrices 



Table 1. Maximum value of fj, 



d=l 


d= 1.1 


d= 1.5 


d = 2 


fast delay 


0.384 


0.367 


0.331 


0.313 


0.289 



0<Pi, Pa, Ps, S >0,Yi, Y-2, T, P, Ra such that 
the following Linear Matrix Inequality (LMI) 



r„ 








<0, 



(15) 



where 

r„ = 





At 

-S - T -T^ 



-hT 
-hR 



+ 



P^ 
+ 



I 
T 





P = 



At, 
-1 

Pi 

P2 P3 



P + 



S 
hR 



Y = [YiY2]. 



(16) 

is feasible. Here J^{p) is given by (10) in case C 
and ^(00) = 7/4 in case B. 

LMI (15) is convex in J^{p) and thus the smaller 
values of J^{p) lead to a less restrictive conditions. 

The time-domain criteria give sufficient conditions 
for the frequency domain Theorem 1. 

Example (Kharitonov & Niculescu, 2003): Con- 
sider the system 



xit) = 
r{t) = 



1 

-1 -2 

+ v{t), 





-1 1 

\v{t)\<f^, m<d. 



x{t) 



X{t-T{t)), 



(17) 

In (Fridman, 2004) the maximum value of /i, for 
which the system is asymptotically stable, was 
found to be = 0.271 for all d > 1. The 
latter result was less restrictive than the one 
by (Kharitonov & Niculescu, 2003). By the time 
domain criterion of (Pridman & Shaked, 2005) for 
d = 1 the corresponding value of /x is greater 
(/U = 0.384), while for d > 1 the result is the same 
{jjL = 0.271). Theorem 2 of the present paper leads 
to a wider stability interval for d > 1 (sec Table 
1). Note that the results by Kao & Rantzer (2005) 
do not improve the existing results by descriptor 
approach and do not treat the case of /i — /x > 0. 



4. CONCLUSIONS 

Linear systems with bounded time- varying delays 

are analyzed under the assumption that the nom- 
inal delay values are not equal to zero. Two cases 
of delay are considered: case B (without any con- 
straints on the delay derivative) and case C (where 
the delay derivative is not greater than d> 1). An 
input-output approach to stability of such systems 
is known to be based on the boimd of the Pa-norm 
of a certain integral operator. In the present paper 



for the first time a tight d — dependent bound 
is derived. The existing bound in ease B is also 
improved. In the past, case C was treated as case 
B, which was restrictive. The new bounds lead 
to improved stability criteria and gives tools for 
further improvements. 



5. APPENDIX 

Proof of Lemma 2. Denote by </? : ^ {0, 1} 
the characteristic function of the domain D in the 
positive quadrant, bounded by the line s = t — h 
and by the graph of the function s = t — h — r]{t), 
i.e., 



1, if {t- h- s){t- h-f]{t) - s) <0, 
0, if {t- h- s){t- h- r]{t) - s) > 

(shaded region in Figure 1). Then z{t) given by 
(8) satisfies the following: 

\\4l,= I i I f{t,si)y{si)dsi 



DC' 







00 



X I / ^{t,s2)y{s2)ds2 ) dt 

nOO pOO / pOO 



(fi{t, Si)ip{t, S2)dt 



'Oo J —00 \^ 

xy^ {si)y{s2)dsids2 



k{si, S2)y'^ {si)y{s2)dsids2, 

k{si,S2)= / f{t,Si)f{t,S2)dt, Si,S2e 



-00 ^ —00 



Hence \\z\\l^ < ||/C||l, • \\y\\l^, where \\JC\\l2 is the 
La-norm of the operator JC : L2{M.) L2{M.) 



k{t,s)f{s)ds, fGL2{ 



By the Riesz-Thorin interpolation theorem (see, 
for example, Okikiolu (1971), Theorem 5.1.3), 
\\IC\\l, < VII^IUi • II^IUoo- Since k{s,,S2) > 
and fc(si,S2) = fc(s2,Si), by the well-known 
formulas for the Pi and Poo-norms, we have 
imiLi = II^IUoc < s^Pse[o,oc)Kis), where 
K{s) = k{si,s)dsi, and hence we decide that 

\\JC\\l, < sup K{s) . (18) 
se[o,oo) 

We shall show that K{t) < 7/4 /i^ for all t € 
[0,00). 

Without loss of generality assume that ri{t) > 0. 
Geometrically, K{t) is the area of the part D{t) of 
the domain D cut out by the strip t—h>S2>t — 
h — r]{t) (double shaded region in Figure 1). 

Thus, D{t) lies inside the parallelogram 

n(t) := {(si, S2) €'M.'^:t-h<S2<t-h- r]{t), 
si — h — < S2 < si — h + fi} 

(sec Figure 2). Divide Il{t) by the vertical lines 
51=^4 and Si = ^5 into two trapezes of total 




p> 1 



t\=t — h — r]{t), t2 =t — h 

Ai : S2 = si - h + /i, A2 : .S2 = .si 
A3 : S2 = Si — /i — /Lt, r : S2 = si - 

Fig. 1. Domains D and D{t) 

So Al 

t2- 



h - ri{si) 




ti h ^5 



ts ^t - II - ri{t), U = t- r]{t), t5 = t, 
te = t + n, t7 = {U + 15)/2 

Fig. 2. Parallelogram Il(t) and domain D{t) 

area /U^ — (/z — and a square (see Figure 

2). The intersection of D{t) with a vertical line 
si = cr, where ^4 < cr < ^5, is contained either in 
the segment [t ~ h — r]{t) , t — h], 01 in the segment 
[t — h,t — h], and hence the area of D{t) Ci {t — 
vit) < si <t} does not exceed 

rt 3 

max{(T - (< - r]{t)), t - rjclr = -v{tf ■ 



L 



t-v(t) 

So we derive the required bound (9) from the 
evident inequality 

which geometrically means that the maximal area 
domain D{t) looks as the shaded region in Figure 
2. 

Proof of Lemma 1. We interpret the function 

T{p) geometrically as follows: 

- for p > 1, T{p)ix^ is the area of the domain 

Dp := {{t, s) e : < s < /z, i - 2/x < s < 
{t- n- s){2s + 2pt - {3p + l)iJ,)) > 0} 

(shaded region in Figure 3(a)), 

- for < p < 1, T{p)iJ? is the area of the domain 

Dp = {{t, s) e R2 : < s < ^, i - 2/x < s < t, 

(t- fx- s){s +pt- 2p/i) < 0} 

(shaded region in Figure 3(b)). 

As in the proof of Lemma 2, we estimate from 
above the value of K{s2)- i.e., the area of the 
domain D{t). 

Without loss of generality we assume that 77 is 
a smooth function, whose zero locus is locally 




p<l 



(a) 




\\ 1 




//////////////// \//A 





(c) 



Fig. 3. The case AT = 1 



finite and consists of only simple roots. Fix t > 0. 

Assume that r](i) > 0, which, in particular, means 
that the line si = t crosses the domain D{t) along 
the segment ti < S2 < ^2 (cf. Figure 2). Put 
N = #{r]-^{0) n{t- II- ri{t),t + 11)) (i.e., the 
number of zeroes in the interval (^3,^5) in Figure 
2). 

Consider few possibilities. 

Step 1. Suppose that A'' = 0. Then D{t) is 
contained in the part of the parallelogram !!(<:), 
right to the line S2 = Si — h (see Figure 2). Then 
its area does not exceed ri{t)n < 11^ < T{p)ii'^. 

Step 2. Suppose that N = 1. The zero r = 
?7~-^(0) n {t — 11 — ri{t),t + 11) corresponds to 
an intersection point of the graphs of the func- 
tions S2 = si — h and S2 = si — h — 
ri{si). This intersection point lies either right 
to the line h := {si = t} or below the line 
h ■■= {S2 = ~psi + {p-l)t- h- ri{t)} (see Fig- 
ure 3(c)). In the former case, the domain D{t) 
remains below the line S2 = si — h, that is we 
have the upper bound from Step 1. In the latter 
case, the part of domain D(t) in the half-plane 

51 > T should lie below the line S2 = si — h and 
above the line I3 = {s2 = — psi + {p — 1)t — h}, 
and the part of D{t) in the half-plane si < r 
should lie below the line ^3 and above the line 

52 = s\ — h (shaded region in Figure 3(c)). It is 
an elementary geometry exercise to show that the 
area of the shaded region in Figure 3(c) does not 
exceed T{p)ii^. 

Step 3. We intend to show that the case N > 1 
reduces to the above considered cases A^ = or 
1. For, we need the following auxiliary geometric 
statement. Consider the parallelogram Il{t), some 
points xi,X2,X3 with decreasing coordinates, ly- 
ing on the line S2 = si — h below the line I2 = 
{s2 = —psi + {p — l)t — h — ri{t)} (see Figure 
4(a)). Draw the lines Z4, Z5 with slope —p through 
the points xi, a;3, respectively, and the vertical line 
le through X2- Denote by F{x2) the area of the 
domain (5(x2), lying inside n(t), between the lines 
^4, ^5, and in the two sectors, bounded by the lines 




Fig. 4. The case N > 1 



Iq and S2 = si — h as marked by dots in Figure 
4(a). Consider F{x2) as a function of the point X2 
running along the segment [a;i, 0:3]. Wc claim that 
F attains its maximum either at X2 = Xi or at 
X2 = 0:3. 

Indeed, when X2 runs along a subscgmcnt of 
[aJijXs] such that the combinatorial type of 5{x2) 
does not change, F is a quadratic polynomial 
in the abscissa t of X2 with the positive second 
derivative. That is, in this subsegment, F attains 
its maximum at an endpoint. If such an endpoint 
differs from Xi and x^, then it corresponds to the 
situation when 

• either X2 belongs to a side of ll{t) (see, for 
example, Figure 4(d)); 

• or the domain ^(^2) intersects with some side 
of n(t) at a point (see. Figure 4(b,c)). 

In the former situation, 5{x2) entirely lies right 
to the vertical line through X2, and then mono- 
tonically grows as X2 tends to a;3. In the latter 
situation, when replacing X2 by xi, the domain 
5{x2) turns into the domain 5{xi) by getting rid 
of 5{x2) n {si < r} and adding the fragment (5_ 
(the trapeze, bounded by the vertical line through 
X2, dashed line, the upper side of n(t), and the line 
S2 = si — h in Figure 4(b)), and when replacing 
X2 by xa, the domain 5{x2) turns into the domain 
5{xz) by by getting rid of 5{x2) n {si > t} and 
adding the fragment 5^ (the trapeze, bounded by 
the vertical line through X2 , dashed line, the lower 
side of n(t), and the line S2 = s\ — h in Figure 
4(b)). One can easy sec that the area of 5(x2) in 
the above situation is less than the maximum of 
the areas of 5{xi) and 5{x3,). The same conclusion 
one can derive in the situation, presented in Figure 
4(c). 



Step 4- Suppose that N > 1. We can take N 
to be odd, adding if necessary one more zero 
close to t — iJ, — r}{t). Denote by xx,X2, ■■■,Xn the 
corresponding intersection points of the graph of 
,S2 = .s'l — h — ri(si) with the line ,S2 = si — h, 
numbered by the decreasing coordinates. Through 
each point X2i-i we draw a line with slope —p, 
l<i<(n + l)/2, and through each point X2i we 
draw a vertical line Z", 1 < i < n/2. Thus, D{t) 
is the union of the regions in 11 (t), bounded by 
the introduced lines as follows (marked by dots in 
Figure 4(e)): 

• above the line l[ and below the line S2 = Si — 
h, 

• below the line l'^, right to the line I'/, and 
above the line S2 = Si — h, 1 < i < n/2, 

• left to the line I", above the line , and 
below the line S2 = si — h, 1 < i < n/2, 

• below the line '(„_|_i)/2 above the line 
S2 = si - h. 

Using the statement of Step 3, we can move the 
point X2 cither to the position xi, or x^ and 
increase the area of S{t). On the other hand, each 
of these limit positions for X2 means that we, in 
fact, have reduced two zeroes of rj in the interval 
{t — fi — r]{t),t + Thus, we inductively come to 
the case N = 1 treated in Step 2. 
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